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Abstract 

Exclusive quarkonium pair production in electron-positron collisions is studied in non- 
relativistic QCD. The obtained cross section for J/tp rj c production in the leading order 
is confronted against the recent measurements by the Belle Collaboration at KEKB. It 
is shown that a large renormalization i^-factor is necessary to explain the experimental 
data. We point out that the J PC = h nature of the hadronic systems that are 
assigned to be r\ c should be tested by the triple angular distributions in terms of the 
scattering angle, and, polar and azimuthal angles of J /tp into leptons. We further study 
J /tp J/ip and T T productions at LEP energies. Although the axial- vector couplings 
of the Z-boson to charm and bottom quarks allow production of such pairs when one of 
them is polarised transversally and the other longitudinally, we find that the integrated 
luminosity at Z pole accumulated by LEP is not large enough to observe the exclusive 
pair production of quarkonium. 
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1 Introduction 



Quarkonium production and decays have long been considered as an ideal means to inves- 
tigate the bound-state formation in QCD. Due to its approximately non-relativistic nature, 
the description of the heavy quark and anti-quark system is one of the simplest applications 
of QCD. For example, the calculation of quarkonium leptonic decays render experimental 
results with a high precision, which may play a crucial role in investigating various phenom- 
ena such as measuring the parton distribution, detecting the Quark-Gluon-Plasma signal 
and even new physics. 

While the quarkonium physics has been studied for more than twenty years, the recent 
interest in the field has been focused on the colour-octet scenario jT] which was triggered by 
the high-pr J/ip surplus production discovered by the CDF collaboration at the Tevatron 
in 1992 [21 El HI- It was proposed based on a novel effective theory, the non-relativistic 
QCD(NRQCD) t 5|. Having achieved the first-step of explaining the CDF data, the colour- 
octet mechanism(COM) had a strong impact into the quarkonium physics and various efforts 
have been made to confirm this mechanism. Although the theoretical framework seems to 
show qualitative agreements with experimental data, there are certain difficulties in the 
quantitative estimate of the colour-octet contribution in particular, in HERA physics 
[7j. It was in such circumstances that the B factory experiments reported their first result 
on the prompt charmonium production from e + e~ collider at y/s = 10.6 GeV jEHHl- As far 
as hadronic uncertainty is concerned, the B factories would provide clearer information of 
the quarkonium production. 

The first result for the inclusive e + +e~ — > J/if>+X process from Belle (with 32.7/6 -1 data 
set) indicated a discrepancy from the theoretical prediction 0. The e + + e~ — > J /if) + cc 
process seemed to dominate the threshold region (z — > 1) of the energy spectra of the 
differential cross section comparing to the COM process e + + e~ — > J/ip+g and colour-singlet 
process e + + e~ — > J /if) + gg, contrary to the theoretical expectation in [101 HU 112 EH Ell- 
in their second report (based on 41.8/fe -1 data set) ^Hj, the direct measurement of the 
e + + — ► J/ip + cc process is presented simultaneously and it is found that experimental 
result is about 10 times larger than the theoretical prediction for this process. In ^3], the 
total cross section of the exclusive e + + e~ — > J /if) + i] c process is found to be: 

a(e + + e~ -> J /if) + rj e ) x B(r] c ^> 4charged) = (0.033±g;gg| ± 0.009)pb (1) 

Note that the observed number of events for J /if) rj c is (67j^]|). Motivated by this measure- 
ment, we investigate J /if) r\ c production at e + e~ colliders. 

The reminder of the paper is organized as follows. In Sec. 2, we will give all the formulae 
used in our analysis. Motivated by the observed large cross section of the J /if) r) c production 
by Belle, we further consider the heavy quarkonium pair production at LEP energies. For 
this purpose, we include the formulae for the Z intermediated processes. In Sec. 3, our 
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numerical results are presented. The double J/ip production through 7*7* intermediated 
states which has been proposed as an explanation of the large cross section of the J /if) r\ c 
process ^H] is also examined. Finally we give our summary and conclusions in Sec. 4. 



2 Formulae 



In this section, we give the formulae which we use in the following sections: the e + e~ 
annihilation into a pair of IS 1 charmonium- and bottomonium-sates. We use the standard 
method in our calculation: we start from the double cc (or bb) production amplitudes, and 
project out the heavy quark and anti-quark pairs into the S-wave states in the colour-singlet 
(see Fig. 1). 

The spin projection operator for the quarkonium production is given by 

P P 11 

Ps,s z { p , 0)= X, v (j - ?! s 2)u{— + q; si) < -, si, -, s 2 \S, S z > , (2) 

where P and S, S z are respectively the quarkonium four-momentum, its spin and the z 
component of the spin; q is the relative momentum of the heavy quarks; and si, s 2 represent 
their spins. In the non-relativistic limit, the covariant forms of the projection operators are 
very simple: 

V ,o(P;0) = ^ 5 (P + M) (3) 

V hSz (P;0) = ^f(P,S z )(P + M) (4) 

respectively, for the pseudoscalar and the vector quarkonium. Here e^(P,S z ) denotes the 
polarization vector of the spin-1 quarkonium state, and M = 2m is the mass of the quarko- 
nium. Projectors (J3J) and (@J) map a QQ pair into the S-wave states. In the rest frame of the 
vector meson, the polarization vector is given by 



e£ = e^(P; S z = 0) = (0, 0, 0, 1) , 4 = ^{P\ S z = ±1) = (0, =F, -i, 0)/V2 , (5) 

for the longitudinal and transverse polarizations, respectively. We need to boost the polar- 
isation vector from rest frame to the laboratory system along the spin quantization axis so 
that S z denotes the helicity (A) in the laboratory frame. In the frame where the electron 
beam is along the z-axis and the quarkonium scattering angle is 9, the polarization vectors 
read 

e £ = 7 (/3, sin 9, 0, cos 9) , 4 = (0, T cos 9, -i, ± sin 0)/V2 . (6) 
respectively, where 

/ 4M 2 1 
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+2 flipped diagrams 

Figure 1: Feynman diagrams for the J /if) rj c production from e + + e~ annihilation. 

Now let us show our analytical formulae for the cross sections. The helicity amplitudes 
for the process e + + — ► J/ip + r] c (T + r] b ) is given by 



_ Q4g 2 e 2 R Q (0Y 

MX=± ~ 37TS3/2 



V Q (v e ~ (TOe) 

s — m| + irrizTz 



^ V el a P VQ pP nQlx f v (8) 

Z "T" ""Z 1 Z J 

M a=o = (9) 

where a denotes the electron helicity in units of 1/2, A the J/ip (T) helicity, and Q is c or 6. 
The initial current is given as j% = (0, —i, a, 0) and Py Q and P^ Q are the four momentum of 
vector and pseudoscalar mesons, respectively. The vector and axial-vector couplings of the 
Z-boson are 

l-4sin 2 ^ 1 J 3 QL - 2e Q sin 2 6 W J 3 Qi 



4 sin 9 W cos 4 sin 0^ cos 6*^ 2 sin cos 6* w 2 sin 6*^ cos 

(10) 

where the is the third component of the weak isospin of the left-handed quark doublet; 
eg is the Q charge in units of proton charge. The absence of the A = amplitudes, the 
scattering angle dependence of the A = ± amplitudes as well as their relative sign tell us the 
spin parity (J p = 0~) of the hadronic system assigned as rj c . 

The above predictions can be tested by experiments through the decay angular distribu- 
tions of Vq —> l + l~, {V c = J/ip, Vb — T). Using a definition D% = M(V\ — > l a d a i) where a' 
is the l~ helicity in the m\ = limit (I = e, /i), we obtain 



Dt = \ —(a' ± cos9*)-^e Tia '^ (11) 
± V167T V ' \/2 V ; 



for A = ± and 



< = fe sinr (12) 

for A = 0. The normalization for the decay amplitudes Eqs. (jllj) and fll2j) is 

2 y dcos0*<ty*|Df | 2 = £ = B{V Q -> fT) (13) 
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Here 9* and 0* are the polar and azimuthal angles of l~ in the Vq rest frame. The 9* 
is measured from the Vq momentum direction in the e + e~ collision rest frame, and 0* is 
measured from the scattering plane ((f)* = tt/2 is along the k x q direction where k and q 
are, respectively the electron and Vq three momenta in the e + e~ collision cm. frame). The 
triple angular distributions are then obtained as 



da 



d cos 9d cos 9*c 



era A=± 



64\ 2 vraV 



r p 2 



35 



2e Q v Q v e (s - ra|) 



2 \2 



+ (m z T z ) 



4/33 



l^(o) \y 



16tt 



;i + cos 2 0)(1 + cos 2 0*) - sin 2 9 sin 2 0* cos : 



(14) 



where the initial beams are unpolarized and final lepton polarizations are unobserved. The 
cos 9* dependence tells us that only the transversally polarized Vq are produced and the 0* 
dependence tells us the relative phase of the A = + and A = — amplitudes. These predictions 
of the Vq tjq production processes should be tested experimentally. After integrating out 
the Vq decay angles and the scattering angle, we find the total cross section 



da 



dcos9 



64 
~3~ 



2 na 2 a 2 



e h 2e Q v Q v e (s 



rrir 



SVP 



[al + v 2 ) 



[s - m 2 z ) 2 + (m z r z y 



|i? Q (0)| 4 /3 3 (l + cos 2 i 



(15) 



We note in passing that a pair of Vq (e.g. J/ip + J ftp or T + T) cannot be produced from 
a single 7* state because of charge conjugation invariance. Such pairs can be produced via 
27* intermediate states or in the Z-boson decays through its axial-vector couplings to the 
heavy quarks. The former process has been studied in Ref. [TH]. For the latter process, we 
find the helicity amplitudes 

32g 2 e 2 a Q R Q (0) 2 M (c 



Mix, 



av. 



(pH _ pfJ- \ AC 



(16) 



37rs 3 / 2 (s — m\ + im z T z ) 

where a is the electron helicity (in units of 1/2), Ai and A2 are the helicities of Vq with 
cos 9 > and cos 9 < 0, respectively, j£ = (0,ia, 1,0) is the initial e + e~ current, and Pv Q i,2 
is the momentum of Vq's. From this result, one can easily find that the amplitudes for 
(Ai, A 2 ) = (±, ±), (0, 0), (±, =f) are zero. The absence of the Ai = A 2 = ± amplitudes is in 
accordance with the Yang's theorem (that forbids the transition between the spin 1 state and 
a pair of identical massless vector bosons), the Ai = A 2 = amplitude is absent due to Bose 
symmetry 17] and the Ai = — A 2 = ± amplitudes vanish because of angular momentum 
mismatch. It is only a pair of longitudinally and transversally polarized Vq's that can be 
produced via the Z-boson exchange in e + e~ annihilation. 



The differential cross section after summing over the Vq helicities for unpolarized beams 



is 



da 



dcos9 



32\ 



2 na 2 a 2 



aQ( a l + v e) 



s 2 (s — m 2 z ) 2 + (m z T z ) 2 



|i? Q (0)| 4 P 5 (1 + cos 2 



(17) 
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From Eqs. (JBJ) and ()16|) . we find the Z-boson decay widths: 



r( z^) = (D' j^yrg (18) 



z 

4/35 



32^ 2 2aa 2 a 2 3 | J R Q (0)| 4 /? 

5~ 

Z 



T(Z - W = ( 3-j %f (19) 



3 Numerical Results 

The following input parameters are used in our numerical analysis in this section: 

a = T | 7 , I c =\, I b = -|, M z = 91.2GeV, T z = 2.495GeV, sin 2 9 W = 0.231, 

Tj,^ e+e - = (5.26 ± 0.37) x 10~ 6 GeV, r T ^ e+e - = (1.32 ± 0.05) x 10~ 6 GeV (20) 

M JH> = 3.1GeV, M T = 9.5GeV 

Although our analysis is strictly in the leading order of perturbative QCD, we adapt the 
two-loop running coupling constant of MS scheme in order to define the "K" factor between 
the leading order prediction and the observed cross section. More specifically, we adopt 

1 ffilnln(/j 2 /Ay , . 

4tt Wms /?oln(/i 2 /Ay /? 3 ln 2 (/i 2 /Ay 1 1 

where ffi = (33 - 2iV / )/3,/?f / = 102 - 10JV) - 8JV//3. a s {M z ) m = 0.118 corresponds to 
Ams — 0.226 GeV for Nf = 5. We set \i = y/s/A as our leading order estimate, because 
the invariant mass of the exchanged gluons in the diagram of Fig. 1 is \i — y/s/2, and 
we account for the factor of two mismatch between the momentum scale and the MS scale 
[TH] . The value of the wave function at the origin can be extracted from the leptonic widths 

r(^-z+r): 

1^(0) | 2 = (22) 

\R b (o) I 2 = ^JrWe- (23) 

Using the experimental values [19J, we obtain: 

|i? c (0)| 2 = (0.53 ± 0.04)GeV 3 , \R b (0)\ 2 = (5.0 ± 0.2)GeV 3 . (24) 

First of all, we show our numerical result for the total cross section of e + + e~ — ► J/ip + 1] C 
which has already been measured by the Belle collaboration. Using the central values of Eq. 
and Ams = 0.226 GeV, Eq. gives 

a(e + + e" -> J /if) + rj c ) = 0.0023pb. (25) 
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Figure 2: The yfs dependence of the total cross section for the exclusive IS double-quarkonium 
production (in Log-scale) integrated over | cos#| < 0.9 is shown. The i^-factor defined in Eq. I|30|) 
is multiplied for the J /if) rj c and T ijb production processes and the J /if J /if and T T processes 
from the intermediated Z boson. 



If all the hadronic systems are indeed rj c decay products and if we set 

B(Vc ^> 4charged hadrons) = 0.6 ± 0.1 (26) 

we find: 

K{VS = 10.6GeV) = g (f + e -^^ + ^W = (24 ± 12) (2?) 

where the error includes only experimental ones, Eq. $IJ) and that in our estimate Eq. (|26|1. 

Let us examine if the large i^-factor in Eq. (|2T|) can be explained. It is well known that 
the leptonic width formulae Eqs. (}2*2"|) and ()23|) suffers from large NLO corrections. As a 
result, our estimate of the value of the wave function at the origin could include large error. 
For instance, one of the potential model calculations give [201 12*T] 

\R C (0)\ 2 = 0.810GeV 3 , |i4(0)| 2 = 6.477GeV 3 . (28) 

By replacing flMJ) by (H, we have a factor of (81/53) 2 . If we use Ams = 0.296 GeV for 
as{M z ) = 0.123 and change the scale fi = y/s/8 instead of our standard choice of [i = a/s/4, 
we find as{^ = y/s/8) = 0.41 instead of as(/i = y/s/A) = 0.29. This gives rise to another 
factor of (35/26) 2 . The product of the two factors is about 5.8 which is still significantly 
smaller than the the value Eq. Q indicated from the experiment. There should be further 
large contributions in the production amplitude and/or the hadronic system should contain 
significant amount of non-7/ c contributions. 

In Ref. [HI], it has been shown that the J /if) pair production process via 7* 7* in- 
termediate state (e + + e~ — > 7*7* — ► J /if) + J/tf>) should be as large as (or larger than) 
e + + e _ — > J /if; + rj c . Because J /if) — > hadrons can not be distinguished easily from 
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rjc — > hadrons, the Belle data Eq. ([IJ may contain contributions from J /if) pair produc- 
tion. The cross section is estimated as: 

a(e + + e~ -> J /if) + J /if)) = 0.0027pb (29) 

Only a small fraction of the observed cross section of Eq. ((TJ can come from the J / if) pair 
production process, because the normalization of the prediction Eq. (}29|) has little theoretical 
uncertainty. We note here that the J / if) pair production process has the following distribution 
properties: the pair of J /if) are transversally (either (+, — ) or (—,+)) polarized and the 
differential cross section behaves as (1 + cos 2 #)/(l — cos 2 9) while the J /if) r\ c production 
behaves as (1 + cos 2 0) (see Eq. (|T5j)). The absence of interference between A = + and 
A = — amplitudes predicts that there is no azimuthal angle (</>*) dependence. The overall 
normalization of the J/ if) pair production contribution should soon be found experimentally 
once the double leptonic decays of the J/ if) pair are observed. 

Our numerical results are summarized in Fig. 2. In this plot, we multiplied the cross 
section for J /if) r\ c production by the K-factor defined: 

2 

(30) 

We use the same K factor, with i^(10.6 GeV) = 24, for the single 7* and Z* exchange 
contributions to J /if) J /if), T 775, and T T production processes, even though the K factor 
for the bottomonium pair production may be smaller than that for the charmonium pair 
production. We hoped that we could study J/if> r\ c production and related processes at Z 
pole by using the LEP data. However considering the LEP integrated luminosity at Z pole 
(about 1 fbr 1 ) and the branching ratio of the leptonic decay of J/if> (0.06), it is unfortunately 
impossible to observe this process at LEP. We find for the input parameters of Eq. (ffilj) and 
K(Mz) = 9.47 from Eq. (J30|) the following branching fractions: 

3.39 x 10~ 13 (31) 

5.73 x 10^ 13 (32) 

9.23 x 10~ n (33) 

4.61 x 10" 11 . (34) 

Even with the cut-off of | cos 6*| < 0.9, the J /if) pair production from two virtual photons 
dominate over all the other exclusive charmonium and bottomonium pair production pro- 
cesses at all energies except around the B factory energies. This is essentially because the 
form factor of the exclusive heavy quarkonium production process drops sharply as s~ 3 / 2 at 
high energies, as can be seen from Eq. (JT3J): 



q a (ys/4) 

"s(\Ao/4) 



B(Z^J/^ + Vc ) = 

B(Z J/1P + J/1P) = 

B(Z^T + Vb ) = 

B(Z -^T + T) = 



An extra suppression factor of s _1 as compared to the high-energy behavior of the light meson 
pair production processes reflects the non-relativistic constraint that the two constituents 
should have the same velocity. 

4 Conclusions 

Exclusive J /if) 7] c production in e + e~ collisions at y's = 10.6 GeV is studied in view of the 
recent Belle observation [Ej- The observed total cross section turns out to be more than one 
order of magnitude larger than the naive leading order prediction of non-relativistic QCD, 
resulting in the huge renormalization factor of K = 24 ± 12. We find that the K factor of up 
to about 6 can be obtained by taking account of the next-leading-order corrections to the 
J /if} leptonic width and by studying the scale dependence of the leading-order prediction. 
On the other hand, the experimental data may contain contributions from non-?7 c origin 
hadronic events, such as hadrons from J /if) decays in the J/if)-pa,ir production via two virtual 
photon exchange JE], and hadrons from two gluon jets in the color-singlet J /if) + gg process 
[TIH EH U21 EH HI] • We propose to use the triple angular distribution of the J/ if> production 
and J I if) decay into charged leptons to test the J/tf> + r) c hypothesis. A peculiar azimuthal 
angle dependence of the lepton distribution about the scattering plane is predicted. 

We have also studied J /if) + J /if) production via Z-boson exchange and find that a pair 
of a transversally polarized J /if} and a longitudinally polarized J /if) can be produced in Z 
boson decays via its axial- vector coupling to the charm quark. Unfortunately, the branching 
fraction of the Z boson decays into J/tf> + J /if}, J/ip + r] c , T + T, T + rjb are all too small to 
be observed in the LEP data, even with a possible large K factor. 

Before closing this report, we note that pair production of S"-wave and P-wave charmo- 
nium has been studied in Ref. [221 I2H| and additionally S'-wave + P-wave as well as P-wave 
+ P-wave charmonium productions have been studied in Ref. [22]. We confirm their results 
of J /if) Xcj (J=0,l,2) production cross sections. Although the perturbative calculation of 
the J I if) r] c production cross section falls short of the observed one, it is still interesting to 
test whether the ratio among cross sections of all the above processes are consistent with 
the predictions of NRQCD. 
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